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DESCRIPTION 

T^/rc^rg^^ rAT.rm.ATTNG A WIRE PA CKING, APPARATUS THEREFOR , AND 

PROGRAM THEREOF 

Technical Field 

The present invention relates to a method of and an 
apparatus for calculation for obtaining an outside diameter of 
a wire packing by bundling and packing a plurality of wires into 
the smallest possible circular shape by using a computer. 

Background Art 

A wire-like structure called a wire harness, which is 
formed by bundling a plurality of wires such as electrical wires 
for electrically connecting electronic equipment, electronic 
components, and the like, is routed in a vehicle or a building 
structure. In recent years, from the viewpoint of improvement 
of space efficiency and the like, such a wire harness has come 
to be required to be as compact as possible without causing a 
decline in electric characteristics. In consequence, there is 
a need to calculate the outside diameter of the wire harness 
more accurately. Conventionally, however, that calculation is 
conventionally performed by experience or a method of 
calculation shown below, and particularly effective 
calculation method has not been proposed. 
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Referring to Figs . 8A and 8B, a problem in the conventional 
method of calculating the outside diameter of a wire harness 
will be described. Fig. 8A shows a wire harness which is made 
up of a relatively small number of wires, and Fig. 8B shows a 
wire harness which is made up of a relatively large number of 
wires . 

In the conventional method of calculation, as shown in 
Fig. 8A, when n wires ai are given, the cross section of these 
wires ai, or a wire harness 10 consisting of these wires ai, 
is assumed to be circular. A diameter L of the wire harness 
10 consisting of these wires ai is determined on the basis of 
the formula: 

7cL 2 /4 = K2(7ili 2 /4) 

Here, 1± denotes the diameter of the wire ai, and K denotes 
a gap coefficient. Since it is difficult to determine 
respective gap coefficients in such a manner as to fit the number 
and layout of various wires, a fixed value is always used for 
all the cases irrespective of the number and layout of wires 
making up the wire harness 10. 

However, as shown in Figs. 8A and 8B, even in the case 
of the wire harnesses having equal diameters L, it is evident 
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that areas of their gaps 11 and 11' differ depending on the number 
and layout of the wires ai and a'i making up the respective wire 
harnesses. If that diameter L is nevertheless determined by 
always using the gap coefficient K of a fixed value in the 
above-described manner, that value becomes inaccurate, so that 
this method has been impractical* 

Thus, there has been no practical method for determining 
the outside diameter of a wire harness which is made up of an 
arbitrary number of wires. It goes without saying that there 
has been no effective method of calculation for obtaining the 
outside diameter of the wire harness when the wires are bundled 
and packed into the smallest circular form possible. 

Disclosure of the Invention 

Accordingly, in view of the above-described 
circumstances, an object of the invention is to provide an 
effective method of calculation for obtaining an outside 
diameter of a wire packing by bundling and packing an arbitrary 
number of wires into the smallest possible circular shape, an 
apparatus and a program therefor. 

To attain the above object, in accordance with a first 
aspect of the invention there is provided a method of 
calculating a wire packing for obtaining an outside diameter 
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thereof by bundling and packing a plurality of wires into the 
smallest possible circular shape so as not to overlap each other 
by using a computer, characterized by comprising: an 
including-circle assuming step of assuming an including circle 
which includes a plurality of circles arranged in a plane so 
as not to overlap each other by assuming that cross-sectional 
shapes of the plurality of wires are the plurality of circles 
having diameters corresponding to respective outer shapes 
thereof; a target-circle defining step of determining a target 
circle which has the same center as that of the including circle 
and is slightly smaller than the including circle, and from 
which at least one of the plurality of circles protrudes; a 
searching step in which the circle protruding from the target 
circle is set as an insertion trial circle, and positions are 
searched to which the plurality of circles other than the 
insertion trial circle can be moved as distantly as possible 
within the target circle without overlapping each other; an 
inserting step of inserting the insertion trial circle in a 
space in the target circle created by changing the layout of 
the plurality of circles on the basis of a result of search in 
the searching step; and a first search controlling step in which 
in a case where all the insertion trial circles have been 
inserted in the target circle, a new target circle which is 
slightly smaller than a present one and has the insertion trial 
circle is set, and the operation returns to the searching step, 
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wherein the including circle is made gradually small by 
repeatedly executing the target-circle defining step, the 
searching step, the inserting step, and the first search 
controlling step. 

In accordance with the above-described first aspect of 
the invention, the including circle surrounding the plurality 
of wires is made gradually small by repeatedly executing the 
target-circle defining step, the searching step, the inserting 
step, and the first search controlling step by using a computer. 
Consequently, the outside diameter of the circle surrounding 
the plurality of wires is obtained efficiently. 

The method of calculating a wire packing in accordance 
with a second aspect of the invention devised to attain the above 
object is characterized in that, in the method of calculating 
a wire packing according to the first aspect, in the searching 
step, a circle Voronoi diagram is constructed by a circle set 
excluding the insertion trial circle and one of the plurality 
of circle, and the target circle, and an examination is made 
with respect to the plurality of circles other than the 
insertion trial circle as to whether or not a center of the one 
circle tangent to both side circles forming each of boundary 
edges in the circle Voronoi diagram is located on the boundary 
edge, to thereby search positions to which the circles can be 
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moved within the target circle. 

In accordance with the above-described second aspect of 
the invention, the search of movement candidate positions of 
the insertion trial circle is extremely simplified by making 
use of a circle Voronoi diagram. 

In the method of calculating a wire packing in accordance 
with a third aspect of the invention devised to attain the above 
object, the method of calculating a wire packing according to 
the first or second aspect is characterized by further 
comprising: a second search controlling step which is 
repeatedly executed together with the target-circle defining 
step, the searching step, the inserting step, and the first 
search controlling step, and in which in a case where the 
insertion of the insertion trial circle is impossible, the 
operation returns to the searching step after determining a new 
target circle which is of an intermediate size between the 
including circle and the present target circle and which has 
the insertion trial circle. 

In accordance with the above-described third aspect of 
the invention, a second search controlling step is further 
comprised in which in a case where the insertion of the insertion 
trial circle is impossible, the operation returns to the 

6 
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searching step after determining a new target circle which is 
of an intermediate size between the including circle and the 
present target circle. Consequently, the outside diameter of 
the circle surrounding the plurality of wires is obtained more 
5 efficiently. 

To attain the above object, in accordance with a fourth 
aspect of the invention there is provided an apparatus for 
calculating a wire packing for obtaining an outside diameter 

10 thereof by bundling and packing a plurality of wires into the 
smallest possible circular shape so as not to overlap each other, 
characterized by comprising: including-circle assuming means 
for assuming an including circle which includes a plurality of 
circles arranged in a plane so as not to overlap each other by 

15 assuming that cross-sectional shapes of the plurality of wires 
are the plurality of circles having diameters corresponding to 
respective outer shapes thereof; target-circle defining means 
for determining a target circle which has the same center as 
that of the including circle and is slightly smaller than the 

20 including circle, and from which at least one of the plurality 
of circles protrudes; searching means in which the circle 
protruding from the target circle is set as an insertion trial 
circle, and positions are searched to which the plurality of 
circles other than the insertion trial circle can be moved as 

25 distantly as possible within the target circle without 



WO 2004/013778 



PCT/JP2003/009821 



overlapping each other; inserting means for inserting the 
insertion trial circle in a space in the target circle created 
by changing the layout of the plurality of circles on the basis 
of a result of search by the searching means; first search 
controlling means in which in a case where all the insertion 
trial circles have been inserted in the target circle, a new 
target circle which is slightly smaller than a present one and 
has the insertion trial circle is set, and the search by the 
searching means is then effected; input means for inputting 
initial information on the plurality of wires; and output means 
for outputting at least the outside diameter of the including 
circle . 

In accordance with the above-described fourth aspect of 
the invention, initial information on the plurality of wires 
is inputted by the input means, the including circle surrounding 
the plurality of wires is made gradually small by the 
target-circle defining means, the searching means, the 
inserting means, and the first search controlling means, and 
the outside diameter of that including circle is outputted by 
the output means. Consequently, the outside diameter of the 
circle surrounding the plurality of wires is obtained 
efficiently. 



The apparatus for calculating a wire packing in 
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accordance with a fifth aspect of the invention devised to 
attain the above object is characterized in that, in the 
apparatus for calculating a wire packing according to the fourth 
aspect, the output means outputs position information on the 
including circle and the plurality of circles. 

According to the above-described fifth aspect of the 
invention, since the position information on the including 
circle and the plurality of circles is outputted, the outside 
diameter of the circle surrounding them, including the layout 
of the plurality of wires, is obtained efficiently. 

The apparatus for calculating a wire packing in 
accordance with a sixth aspect of the invention devised to 
attain the above object is characterized in that, in the 
apparatus for calculating a wire packing according to the fourth 
or fifth aspect, the searching means includes second searching 
means in which a circle Voronoi diagram is constructed by a 
circle set excluding the insertion trial circle and one of the 
plurality of circle, and the target circle, and in which an 
examination is made with respect to the plurality of circles 
other than the insertion trial circle as to whether or not a 
center of the one circle tangent to both side circles forming 
each of boundary edges in the circle Voronoi diagram is located 
on the boundary edge, to thereby search positions to which the 
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circles can be moved within the target circle. 

In accordance with the above-described sixth aspect of 
the invention, the search of movement candidate positions of 
5 the insertion trial circle is extremely simplified by making 
use of a circle Voronoi diagram in the second searching means. 

In the apparatus for calculating a wire packing in 
accordance with a seventh aspect of the invention devised to 
attain the above object , the apparatus for calculating a wire 
packing according to any one of fourth, fifth, and sixth aspects 
is characterized by further comprising: second search 
controlling means which, in a case where the insertion of the 
insertion trial circle is impossible, causes the search by the 
searching means to be effected after determining a new target 
circle which is of an intermediate size between the including 
circle and the present target circle and which has the insertion 
trial circle. 

20 In accordance with the above-described seventh aspect of 

the invention, in a case where the insertion of the insertion 
trial circle by the second search controlling means is 
impossible, the search by the searching means is effected again 
after determining a target circle which is of an intermediate 

25 size between the including circle and the present target circle . 

10 
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Therefore, the outside diameter of the circle surrounding the 
plurality of wires is obtained more efficiently. 

Brief Description of the Drawings 
5 Fig. 1 is a block diagram illustrating an example of a 

hardware configuration for realizing a calculation method and 
an apparatus therefor in accordance with the invention; 

Fig. 2 is a flowchart illustrating a basic processing 
procedure in accordance with an embodiment of the calculation 
10 method of the invention; 

Fig. 3A is a diagram illustrating an initial state; 
Fig. 3B is a diagram illustrating an insertion trial 
circle protruding from a target circle; 

Fig. 3C is a diagram illustrating a state in which the 
15 insertion trial circle shown in Fig. 3B is inserted in the target 
circle; 

Fig. 3D is a diagram illustrating a final result; 
Fig. 4 is a flowchart illustrating search and insertion 
processing shown in Fig. 2; 
20 Fig. 5 is a flowchart illustrating first search 

processing; 

Fig. 6 is a flowchart illustrating second search 
processing; 

Fig. 7A is a diagram illustrating an example of a circle 

25 set ; 

11 
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Figs . 7B and 7C are a circle Voronoi diagram and a Laguerre 
circle Voronoi diagram, respectively, with respect to the 
circle set shown in Fig. 7A; 

Fig. 8A is a diagram illustrating a wire harness which 
5 is made up of a relatively small number of wires; and 

Fig. 8B is a diagram illustrating a wire harness which 
is made up of a relatively large number of wires. 

Best Mode for Carrying out the Invention 

Referring now to the accompanying drawings, a description 
will be given of an embodiment of the invention. 

First, referring to Fig. 1, a description will be given 
of a hardware configuration for realizing a method of 
calculating a wire packing. Fig. 1 is a block diagram 
illustrating an example of a hardware configuration for 
realizing a calculation method and an apparatus therefor in 
accordance with the invention. 

As shown in Fig. 1, the hardware configuration is realized 
20 by a known personal computer, a general-purpose computer, or 
the like. This computer is comprised of an input unit 1, an 
input/output interface circuit (I/O) 2, a central processing 
unit (CPU) 3, a memory 4, an output unit 5 and a read/ write unit 
6. The input unit 1, the memory 4, the output unit 5 and the 
25 read/write unit 6 are electrically connected to the CPU 3 

12 
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through the I/O 2 and the like. 

The input unit 1 is, for example, a keyboard and a mouse 
device which are used for inputting input data in processing 
5 which will be described later. The CPU 3 includes a control 
section 31 for controlling the input unit 1, the output unit 
5, and the like, as well as an operation part 32 for effecting 
processing concerning the calculation method, which will be 
described later, in accordance with programs stored in the 
10 memory 4. 

The memory 4 includes a program memory 41 for storing 
programs and the like corresponding to each processing in 
accordance with the calculation method which will be described 
15 later, as well as an operation memory 42 to which work areas 
for various processing to be performed by the CPU 3 are allotted. 
The output unit 5 is, for example, a monitor display and a printer 
for outputting the results of processing performed by the CPU. 

2 0 The read/write unit 6 is a device which reads out a wire 

packing calculating program 7a according to the present 
invention which stored in a recording medium 7 such as CD-ROM 
(for example, a program for processing procedures shown in Figs. 
2,4,5 and 6 described later) , and transmits it to the program 

25 memory 41. The read/write unit 6 has a function of writing the 

13 
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calculation result in the recording medium 7 . Incidentally/ the 
computer may includes a communication interface such as 
not- shown modem board and LAN card. 

The CPU 3 installs the wire packing calculating program 
7a read out by the read/write unit 6 in a program memory 41 of 
the memory 4. After power the computer on, the program 7a is 
executed and the computer functions as the apparatus for 
calculating the wire packing. The wire packing calculating 
program 7a is capable of being installed in, for example, a 
personal computer and general-purpose computers having the 
above construction, and when the program is installed, the 
computer functions as the apparatus for calculating the wire 
packing. 

The wire packing calculating program 7a may be provided 
not only through the recording medium such as CD-ROM but also 
through a communication lint such as Internet, dedicated line 
and LAN, for example. 

Next, referring to Figs. 2 to 7, a description will be 
given of a processing procedure concerning an embodiment of the 
calculation method in accordance with the invention. 
Accordingly, referring to Fig. 2, a description will be first 
given of a basic processing procedure of this calculation method 
Fig. 2 is a flowchart illustrating a basic processing procedure 
in accordance with an embodiment of the calculation method of 
the invention. 

14 
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In this calculation method/ the problem resolves itself 
into one in which cross-sectional shapes of a plurality of wires 
making up the wire harness are assumed to be a plurality of 
5 circles having diameters corresponding to the respective outer 
shapes thereof, and when n circular cylinders having these 
circles in their cross sections are bundled, the size of a circle 
capable of surrounding and accommodating the overall circles 
is examined. In practice, an effective calculation method is 
10 conceived for bundling and packing a plurality of wires into 
the smallest possible circular shape and for obtaining its 
outside diameter. 

In the basic processing shown in Fig. 2, radii n, r 2 , . . • , 
15 r n of n circles Ci, c 2 , . . c n whose outer shapes respectively 
correspond to cross-sectional shapes of a plurality of wires 
such as electric wires making up the wire harness, a number p 
which is smaller than and sufficiently close to 1, e.g., p = 
0.95, a terminating reference value 8 which is a sufficiently 
20 small integer, e.g., the terminating reference value s = min ( (ri, 
r2f • r 3 )/100) are given as input information. 

In addition, a radius R of the smallest possible circle 
for packing the n circles Ci, c 2 , . • . , c n so as not to overlap 
25 each other, as well as position information on a circle C and 

15 
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the circles ci, c 2 , . . c n at this time are given as the output 
information. 

For this reason, in Step SI shown in Fig. 2, the circles 
5 Ci, c 2 , . . . / c n are laid out in a plane so as not to overlap each 
other, and a large circle surrounding them, i.e., an including 
circle C, is found. 

Next, in Steps S2, S3, and S4, a circle which has the same 
10 center as that of the aforementioned including circle C and 
whose radius is p-fold that of the including circle C, i.e., 
a target circle D, is determined. Namely, in a loop consisting 
of Steps S2, S3, and S4, a target circle D is determined which 
has the same center as that of the including circle C and is 
15 slightly smaller than the including circle C, and from which 
at least one of the plurality of circles Ci, c 2 , . • . , c n protrudes . 
It should be noted that, in the following processing steps, the 
layout is changed so that the circles ci, c 2 , c n are 

accommodated in the target circle D. 

20 

Next, in Step S5, search and insertion processing is 
effected. Namely, here, circles other than an arbitrary circle 
Ci protruding from the target circle D are removed in order 
starting with one whose distance from the circle Ci is largest. 
25 Those circles which can be placed at farther positions are moved 

16 
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as distantly as possible, and when such movement is impossible, 
the circles are left at the present positions . Then, an attempt 
is made to move the single circle Ci to a space created by such 
movement, i.e., an insertion is attempted. It should be noted 
5 that the processing of this Step S3 will be described later with 
reference to Figs. 4 to 7 . 

Next, in Step S6, a determination is made as to whether 
or not the insertion of the circle Ci has been successful in 

10 Step S5 above. If successful (Y in Step S6) , the operation 
returns to Step S3, and if not (N in Step S6) , the operation 
proceeds to Step S7. When the operation returns to Step S3, 
a determination is made as to whether or not there is any other 
protruding circle. If there is, the search and insertion 

15 processing in Step S5 is effected again with respect to this 
protruding circle. If there is not, the operation proceeds to 
Step S4 to effect the same processing as described above. 

Meanwhile, in Step S7, a circle of an intermediate size 
20 between the aforementioned including circle C and the 

aforementioned target circle D for which the insertion was 
unsuccessful is newly set as the target circle D. Next, in Step 
S8, a determination is made as to whether the difference between 
the radii of the including circle and the target circle D used 
25 in the processing in Step S7 above is equal to or less than the 
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aforementioned terminating reference value e. If this 
difference is greater than the terminating reference value s 
(N in Step S8) , the operation returns to Step S3 to repeat the 
same processing as described above, whereas if this difference 
is equal to or less than the terminating reference value 8 (Y 
in Step S8), the operation proceeds to Step S9. 

In Step S9, the radius of this including circle C is 
outputted to the output unit 5 as the final radius R of the wire 
harness. In addition, the position information on the circle 
C and the circles Ci, c 2 , . . . / c n at this time are also outputted 
to the output unit 5. These outputs may be given on a monitor 
display or may be printed out on paper. It should be noted that 
the aforementioned number p and terminating reference value s 
may be slightly changed, as required. 

The behavior of the respective circles according to the 
above-described processing procedure will be shown with 
reference to Figs. 3A to 3D. Figs. 3A to 3D are diagrams 
illustrating the behavior according to the processing procedure 
shown in Fig. 2. In particular, Fig. 3A shows an initial state; 
Fig. 3B shows an insertion trial circle protruding from the 
target circle; Fig. 3C shows a state in which the insertion trial 
circle shown in Fig. 3B is inserted in the target circle; and 
Fig. 3D is a diagram illustrating a final result. 

18 
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Fig, 3A shows the initial layout of the n circles Ci given 
as well as the including circle C surrounding them. Fig. 3B 
shows a state in which processing is underway, and illustrates 
the target circle D which is slightly smaller than the presently 
obtained including circle C as well as the insertion trial 
circle c n which is one of the circles protruding from this target 
circle D and whose insertion is tried. 

In addition, Fig. 3C shows a state after the search and 
insertion processing in Step S5 shown in Fig. 2 was carried out 
with respect to the insertion trial circle c n shown in Fig. 3B. 
It should be noted that, in Fig. 3C, circles mi (a circle set 
surrounding the coarse hatchings) show a moved circle set which 
was moved to insert the insertion trial circle c n in the 
above-described search and insertion processing. It should be 
noted that, as can be understood from this drawing, there are 
cases where the other protruding circles are also packed in the 
target circle D in the process of the insertion processing of 
the insertion trial circle c n . Then, the result in which 
insertion processing was effected with respect to all the 
protruding circles is shown in Fig. 3D. 

Thus, the layout of the plurality of wires making up the 
wire harness is changed such that they are located as distantly 

19 
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as possible from the wire protruding from the including circle, 
and the protruding wire is inserted in the space thus created. 
As this operation is repeatedly calculated, the outside 
diameter of the wire harness surrounding the plurality of wires 
5 is efficiently obtained. 

Next, referring to Fig. 4, a description will be given 
of the search and insertion processing in Step S5 shown in Fig. 
2 referred to above. Fig. 4 is a flowchart illustrating the 
10 search and insertion processing in Fig. 2. 

In the search and insertion processing shown in Fig. 4, 
radii r± of n circles Ci, their centers (xi, yj , i = 1, 2, 
n, and the target circle D are given as input information. It 
15 is assumed, however, that the n circles c± do not overlap each 
other, and that the final circle c n protrudes from the target 
circle D. There may be other protruding circles. 

Meanwhile, by way of output information, if the final 
20 circle c n can be inserted in the target circle D without causing 
the circles already placed in the target circle D to protrude 
from this target circle D, the central positions of the n circles 
for realizing that insertion are outputted as a successful 
result. If that insertion is impossible, a message to that 
25 effect is outputted as an unsuccessful result. 

20 
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First, in Step S51 in the search and insertion processing, 
the n circles c± are rearranged in the order of distance from 
the aforementioned final circle c n . Specifically, this order 
5 is based on the distance between the center of each of the n 
circles Ci and the center of the final circle c n . Here, the 
circle numbers as the result of the rearrangement are newly set 
as ci, C2, . . . , c n for the sake of simplicity. It should be noted 
that this final circle will be hereafter referred to as an 
10 insertion trial circle. 

Next, processing shown in Steps S52 to S54a (or Step S54b) 
is effected with respect to i = 1, 2, . . . , n - 1 . In Step S52, 
search processing is carried out. Namely, a search is made of 

15 a movement candidate position to which the circle ci can be moved 
without overlapping the other circles within the target circle 
D. Specifically, in this search processing, first search 
processing shown in Fig. 5 or second search processing shown 
in Fig. 6 is carried out. In the first search processing, a 

20 search is made of such a movement candidate position that the 
circle Ci becomes farther from the insertion trial circle c n 
than the present position of the circle c± without overlapping 
the other circles within the target circle D. In addition, in 
the second search processing shown in Fig. 6, by using the 

25 concept of a circle Voronoi diagram, a search is made of a 

21 
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movement candidate position to which the circle c± can be moved 
without overlapping the other circles within the target circle 
D. These searches will be described later. 

5 Then, in Steps S53, S54a, andS54b, if a movement candidate 

position is present in the above-described search processing 
(Y in Step S53) , the circle Ci is moved to a position which is 
farthest from the insertion trial circle c n among such movement 
candidate positions (Step S54a) . If there is no movement 
10 candidate position (N in Step S53) , the circle Ci is left at 
the present position (Step S54b) . After such processing is 
effected with respect to i = 1, 2, . . . , n, the operation proceeds 
to Step S55. It should be noted that the aforementioned Steps 
S52 to S54 correspond to the search process in the claims. 

15 

Next, in Step S55, the insertion of the insertion trial 
circle c n is tried with respect to the space in the target circle 
D created by the loop processing consisting of the 
aforementioned Steps S52 to S54a (or Step 54b) . 

20 

Then in Steps S56, S57a, and S57b, if the insertion is 
successful by the trial of the insertion (Y in Step S56) , the 
insertion trial circle c n is moved to that position (Step S57a) . 
If the insertion is not successful (N in Step S56) , a message 
25 to that effect is outputted (Step S57b) . It should be noted 
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that at the time of the success the central positions of the 
n circles for realizing it are outputted. Then, when the series 
of processing consisting of Steps S51 to S56a (or Step S56b) 
is completed, the operation returns to the ensuing processing 
5 shown in Fig. 2. 

Referring to Figs. 5 and 6, a description will be given 
of two examples of the above-described search processing. 
First, the first search processing will be described with 
10 reference to Fig. 5. Fig. 5 is a flowchart illustrating the 
first search processing. 

In the first search processing shown in Fig. 5, attention 
is paid to the fact that that when there is a space for moving 

15 the n circles ci, in a state in which the circle c ± is moved 
to a most distant position from the aforementioned insertion 
trial circle c n/ the circle Ci should be tangent to two circles. 
However, there are cases where one of these two tangent circles 
is the aforementioned target circle D. Accordingly, here, a 

2 0 set formed by the whole of the given n circles c± and the target 
circle D is set as S = {n circles Ci, c 2 , c n , D} . Then, 

processing shown in the following Steps S521 to S529 is effected 
with respect to all the two circles Cj and c k other than the 
circle c± ^ S. 

25 

23 
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First, in Step S521, two positions are searched where the 
circle ci with a radius r± is tangent to both the circle c^ and 
the circle c^. However, it is assumed that if the circle Cj or 
the circle c k is a circle other than the target circle D, the 
circle c± is tangent to that circle from the outside, and if 
the circle c 3 or the circle c k is the target circle D, the circle 
ci is tangent to that circle from the inside. There are only 
two such positions, and the centers in such a case are set as 
(x' ±/ y'i) and (x"i, y"i) , respectively. 

Next, in Step S522, a determination is made as to whether 
or not if the circle c± is moved to one of the aforementioned 
two positions, the circle c± becomes more distant from the 
insertion trial circle c n than the present position. Namely, 
a comparison is made between a distance X' from one center (x'i, 

o 

y' i) to the center of the insertion trial circle c n and a distance 
X from the center of the circle Ci at the present position to 
the center of the insertion trial circle c n . If the distance 
X' is greater than the distance X (Y in Step S522) , the operation 
proceeds to Step S523; if not (N in Step S522), the operation 
proceeds to Step S52 6 which will be described later. 

In Step S523, an examination is made as to whether or not 
the circle c± with the radius ri, when placed at the 
aforementioned center (x'i, y'i)/ overlaps any of all the 
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circles other than the circle d, the circle Cj, the circle c k , 
and the target circle D. The overlap is determined in Step S524 . 
Here, if it is determined that the circle ci does not overlap 
any circle (N in Step S524) , the operation proceeds to Step S525 . 
5 In Step S525, the aforementioned center (x'i, y' ± ) is added as 
one of the movement candidate positions of the circle ci," 
otherwise (Y in Step S524) , the operation proceeds to Step S526. 

Further, the aforementioned center (x'i, y'i) is replaced 
10 by the other center (x" ± , y'\) , and the following processing 
in Steps S526 to S529 is effected in the same way as Steps S522 
to S525 above. In Step S52 6, a determination is made as to 
whether or not if the circle Ci is moved to the other one of 
the aforementioned two positions, the circle c± becomes more 
15 distant from the insertion trial circle c n than the present 
position. Namely, a comparison is made between a distance X" 
from the other center (x"i, y"±) to the center of the insertion 
trial circle c n and the distance X from the center of the circle 
ci at the present position to the center of the insertion trial 
20 circle c n . If the distance X" is greater than the distance X 
(Y in Step S526) , the operation proceeds to Step S527; if not 
(N in Step S526) , the operation proceeds directly to a next one 
in the loop. 

25 In Step S527, an examination is made as to whether or not 
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the circle Ci with the radius ri, when placed at the 
aforementioned center (x"i, y"i) , overlaps any of all the 
circles other than the circle Ci, the circle Cj, the circle c k , 
and the target circle D. The overlap is determined in Step S528 . 
Here, if it is determined that the circle c± does not overlap 
any circle (N in Step S528) , the operation proceeds to Step S52 9 . 
In Step S529, the aforementioned center (x"i, y"±) is added as 
one of the movement candidate positions of the circle Ci; 
otherwise (Y in Step S528) , the operation directly proceeds to 
a next one in the loop. When such processing is effected with 
respect to all the two circles Cj and c k other than the circle 
c±, the operation returns to the immediately following 
processing shown in Fig. 4. 

By using the above-described first search processing, the 
calculation of the outside diameter of the wire harness, which 
is conventionally performed by the method shown in Fig. 8 or 
by experience, can be improved and performed more accurately. 
Accordingly, an aid is provided in the design of the wire harness . 
On the other hand, there arises a problem in that if the first 
search processing is used, the amount of calculation becomes 
enormously large. Namely, in the first search processing shown 
in Fig. 5 referred to above, since the above-described operation 
is effected with respect to the set of the circle Ci, the circle 
Cj, and the circle Ck, the calculation time becomes O (n 3 ) . In 
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the search and insertion processing shown in Fig. 4, since the 
processing is effected with respect to all of I = 1, 2, . 
n, the calculation time becomes O (n 4 ) . Further, in the 
processing shown in Fig, 2, since processing is effected 
5 repeatedly by incorporating the processing shown in Fig. 4, the 
total amount of calculation becomes enormously large. 
Improvement is made in this respect in the second search 
processing shown in Fig. 6 below. 

10 Fig- 6 is a flowchart illustrating the second search 

processing. Fig. 7A is a diagram illustrating an example of 
a circle, set- Figs. 7B and 7C are a circle Voronoi diagram and 
a Laguerre circle Voronoi diagram, respectively, with respect 
to the circle set shown in Fig. 7A. 

15 

First, a basic concept of this second search processing 
will be shown. In this second search processing, the search 
of a movement candidate position of an insertion trial circle 
d is made efficient by using the concept of a known Voronoi 
20 diagram. Namely, in the above-described first search 

processing, positions where the circle ci is tangent to all sets 
of circles Cj and Ck are determined. However, if the concept 
of the Voronoi diagram is utilized, the candidates can be 
limited. 

25 
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When a finite number of circles not overlapping each other 
in a plane are given, the plane can be divided by determining 
to which circle a given circle is closest . This divided diagram 
is called a circle Voronoi diagram. This circle Voronoi diagram 
5 is also shown in document 1 (A. Okabe, B. Boots, K. Sugihara 
and S.N. Choi: Spatial Tessellations — Concepts and 
Applications of Voronoi Diagrams, 2nd Edition. John Wiley and 
Sons, Chichester, 2000) . 

L0 For example, a circle Voronoi diagram with respect to a 

circle set shown in Fig. 7A is shown in Fig. 7B. A point on 
a boundary edge e 3 which is called a Voronoi edge in Fig. 7B 
has a characteristic that the point is equidistanced from two 
circles c k and ci and the other circles are located more distantly. 
15 Accordingly, a circle which is tangent to the two circles c k 
and ci and does not overlap other circles has a center on the 
boundary edge e 3 of the circle Voronoi diagram. Accordingly, 
as for the movement candidate positions of the circle Ci, it 
suffices to search only the sets of two circles c k and c x 
20 sandwiching the boundary edge ej of the circle Voronoi diagram. 
Since the number of the boundary edges ej of the circle Voronoi 
diagram of n circles is proportional to n, the calculation time 
of the sets of the circles c k and Ci to be searched is 0 (n 2 ) 
in the first search processing shown in Fig. 5 referred to above, 
25 whereas the calculation, time of the sets of the circles c k and 
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ci to be searched is O (n) in this case. 

Further, in the first search processing shown in Fig. 5 
referred to above, a search is made of the overlap with all the 
circles Cj and c k other than the circle c± with respect to the 
movement candidate positions of the circle c ± tangent to the 
two circles c k and ci. Here, however, that search is made 
unnecessary. Namely, it is sufficient if a search is made as 
to whether or not the movement candidate positions of the 
circles Ci are located on the boundary edges ej . This is because 
it can be derived from the characteristic of the Voronoi diagram 
that if the circle c± is located on the boundary edge e j# the 
circle Ci does not overlap another circle, and if the circle 
ci is not located on the boundary edge e j# the circle ci overlaps 
another circle. Accordingly, the calculation time of O (n) can 
be reduced to O (1) . 

The processing procedure of the second search processing 
based on such a concept is shown in Fig. 6 referred to below. 
In Step S521' in Fig. 6, a circle Voronoi diagram is constructed 
with respect to the aforementioned circle set S - {ci}, i.e., 
all the circles other than the circle Ci. Here, since the circle 
set S - {c±} consists of n circles, the number of boundary edges 
is also proportional to n. 
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Then, processing in the following Steps S522' to S525' 
is performed with respect to each boundary edge e 3 ( j = 1, 2, . . . , 
n) . 

In Step S522', a circle c' i with a radius r ± tangent to 
the circles c* and ci on both sides of the boundary edge e-, is 
plotted. Here again, it is assumed that if the circle c k or 
the circle ci is a circle other than the target circle D, the 
circle c'i is tangent to that circle from the outside, and if 
the circle c k or the circle ci is the target circle D, the circle 
c'i is tangent to that circle from the inside. Incidentally, 
there are only two such positions. 

Next, if it is determined in Step S523' that the 
above-described tangent circle c'i with the radius r ± is present 
(Y in Step S523' ) , and if it is determined in Step S524' that 
its center is located on the boundary edge ej (Y in Step S524' ) , 
the operation proceeds to Step S525' to add this circle c'i as 
the movement candidate position of the circle c±. Otherwise 
(N in Step S523' and N in Step S524' ) , the operation proceeds 
directly to the next step. If such processing is effected with 
respect to all the boundary edges e^, the operation returns to 
the following processing shown in Fig. 4. 

As is apparent from the foregoing description, it can be 
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appreciated that the search of movement candidate positions of 
the circles c± can be extremely simplified by utilizing the 
concept of the Voronoi diagram. The aforementioned document 
1 shows that a circle Voronoi diagram of n circles can be 
constructed in the calculation time of O (nlogn) . Accordingly, 
the processing of the aforementioned Step S521' can be executed 
in O (nlogn) . Meanwhile, since the number of the boundary edges 
of the circle Voronoi diagram of the n circles is only 
proportional to n, the processing of the aforementioned Steps 
S522' to S524' can be executed in O (n) time. It follows from 
above that the calculation time of the second search processing 
shown in Fig. 6 becomes O (nlogn) . By way of reference, since 
the calculation time of O (n 3 ) is required in the first search 
processing shown in Fig. 5, it can be seen that the efficiency 

improves substantially. Incidentally, if the second search 
processing in Fig. 6 is incorporated into the search and 

insertion processing in Fig. 4, since the processing in Fig. 

6 is executed O (n) times in the search and insertion processing 

in Fig. 4, the calculation time of the search and insertion 

processing in Fig. 4 becomes O (n 2 logn) . 

It should be noted that a simple method of calculating 
the above-described circle Voronoi diagram is a method in which 
a Laguerre Voronoi diagram described in the aforementioned 
document 1 is constructed, and is then converted into a circle 
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Voronoi diagram by the flip operation of a side described in 
document 2 (D.S. Kim and K. Sugihara: Voronoi diagram of a 
circle set from Voronoi diagram of a point set, I. Topology. 
Computer Aided Geometric Design, vol. 18 (2001), pp. 541 - 562) 
and document 3 (D.S. Kim and K. Sugihara: Voronoi diagram of 
a circle set from Voronoi diagram of a point set, II. Geometry. 
Computer Aided Geometric Design, vol. 18 (2001), pp. 563 - 585) . 

As described above, according to this embodiment, the 
outside diameter of the wire harness surrounding the plurality 
of wires can be efficiently obtained by repeatedly calculating 
the operation in which the layout of the plurality of wires 
making up the wire harness is changed such that, by using a 
computer, the wires are moved as distantly as possible from the 
wire protruding from the including circle, and the protruding 
wire is inserted in the space thus created. In particular, by 
adopting the concept of the circular Voronoi diagram, it becomes 
possible to obtain the outside diameter of the wire harness 
extremely simply and in a short time. 

It should be noted that, in the above-described 
embodiment, the radius is outputted to determine the outside 
diameter of the wire harness, but it goes without saying that 
the diameter may be outputted instead. In addition, an 
arrangement may be provided such that not only is the radius 
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of the wire harness outputted in the above-described manner, 
but information indicating the layout of the respective wires 
is outputted. In addition, the values of p and □ are not limited 
to the values shown in the above-described embodiment, and may 
be modified, as required, within the scope of the gist of the 
invention. 

In accordance with the first aspect of the invention, the 
including circle surrounding the plurality of wires is made 
gradually small by repeatedly executing the target-circle 
defining step, the searching step, the inserting step, and the 
first search controlling step by using a computer. 
Consequently, the outside diameter of the circle surrounding 
the plurality of wires can be obtained efficiently. 

In accordance with the second aspect of the invention, 
the search of movement candidate positions of the insertion 
trial circle is extremely simplified by making use of a circle 
Voronoi diagram. Accordingly, the outside diameter of the 
circle surrounding the plurality of wires can be obtained in 
a short time. 

In accordance with the third aspect of the invention, a 
second search controlling step is further comprised in which 
in a case where the insertion of the insertion trial circle is 
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impossible, the operation returns to the searching step after 
determining a new target circle which is of an intermediate size 
between the including circle and the present target circle. 
Consequently, the outside diameter of the circle surrounding 
5 the plurality of wires can be obtained more efficiently. 

In accordance with the fourth aspect of the invention, 
initial information on the plurality of wires is inputted by 
the input means, the including circle surrounding the plurality 

10 of wires is made gradually small by the target-circle defining 
means, the searching means, the inserting means, and the first 
search controlling means, and the outside diameter of that 
including circle is outputted by the output means. 
Consequently, the outside diameter of the circle surrounding 

15 the plurality of wires can be obtained efficiently. 

According to the fifth aspect of the invention, since the 
position inf ormation on the including circle and the plurality 
of circles is outputted, the outside diameter of the circle 
20 surrounding them, including the layout of the plurality of wires, 
can be obtained efficiently. 

In accordance with the sixth aspect of the invention, the 
search of movement candidate positions of the insertion trial 
25 circle is extremely simplified by making use of a circle Voronoi 
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diagram by the second searching means . Accordingly, the outside 
diameter of the circle surrounding the plurality of wires can 
be obtained in a short time. 

5 In accordance with the seventh aspect of the invention, 

in a case where the insertion of the insertion trial circle by 
the second search controlling means is impossible, the search 
by the searching means is effected again after determining a 
target circle which is of an intermediate size between the 
10 including circle and the present target circle. Therefore, the 
outside diameter of the circle surrounding the plurality of 
wires can be obtained more efficiently. 
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